Flux-induced persistent motion of solitons in topological Josephson junctions 
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We propose a Hamiltonian describing the collective motion of a phase soliton (Josephson vortex) 
within a Josephson junction separating two topological superconductors. Using the Hamiltonian we 
derive the equations of motion for the soliton trapped in an annular Josephson junction. We find 
a universal phase accumulated by the soliton when it encircles the junction, that depends on the 
parity of the number of vortices enclosed by the junction. We demonstrate that the presence of this 
phase can be measured through its effect on the junction's resistance. 



Introduction. — Topological order, a remarkable type 
of quantum order that admits no local order param- 
eter [1-5 , was manifested recently in many forms in 
several laboratories worldwide [6HTU]. Most recently, 
ground-breaking experiments in solid-state devices have 
potentially observed signatures of Majorana fermions [8]- 
[T0] . quasi-particles that are their own quasi- holes [TT] . 
An unambiguous discovery of these highly-prized quasi- 
particles may carry far-reaching consequences for basic 
science and technology alike: they entail a novel type of 
non- commutative quantum statistics [T2HT7] that beck- 
ons a paradigm shift in quantum computation [T8| fT9] . 

Some of the most fascinating constructs that have 
emerged out of these studies are the topological Joseph- 
son junctions (TJJs) [20 , 21 , which differ from their non- 
topological counterparts by the presence of a pair of one- 
dimensional counter-propagating Majorana states nucle- 
ated at the junction. The resulting dynamics for the junc- 
tion is consequently governed by a modified sine- Gordon 
Hamiltonian, where the regular bosonic degrees of free- 
dom couple with the low- lying Majorana fermions [22] . In 
particular, properties of Josephson vortices through the 
junction are modified so that each vortex carries a Majo- 
rana zero mode [20, 22, 24 ; 26 . TJJs have been argued to 
be a promising venue for realizing itinerant non-Abelian 
anyons taking the form of Josephson vortices [22 , 24|[25]. 

In this paper we derive an effective quantum Hamil- 
tonian for a Josephson vortex in a TJJ (see Eq. ([l8|)), 
unveiling a remarkable property of these vortices. It is 
shown that the presence of the low-lying Majorana edge 
states that interact with the vortex crucially affects its 
dynamics. As a consequence, for the case of a Joseph- 
son vortex going around an annular Josephson junction 
[27l [28] (see Fig. [I]), the vortex accumulates a universal 
phase. This phase can be exploited to induce a persis- 
tent motion of the vortex around the junction, triggered 
by the nucleation of an additional vortex in the region 
enclosed by the junction (i.e. by changing the magnetic 
flux $ through the central hole). The induced motion 
drives the Josephson junction into its finite resistance 
state [23] , revealing the presence of the phase. For non- 
topological Josephson junctions, an externally induced 
charge Q can drive the vortex into a similar persistent 
motion [27] through the Aharonov-Casher effect [29H3T] ; 



The system acts as an elegant dual of an Aharonov-Bohm 
ring for electrons. In contrast, for TJJs the persistent 
motion of the Josephson vortex can be controlled through 
two knobs instead of one: i) Continuously, using the in- 
duced charge Q in the region enclosed by the junction; 
and (ii) Using the enclosed flux, whose value is quantized 
in units of h/2e. In units of the charge, the nucleation of 
an extra vortex within the central region is equivalent to 
an e/4 shift in the Aharonov-Casher charge Q (where e 
is the electronic charge). 

We proceed to derive the equations governing the 
motion of a Josephson vortex. Our starting point is 
the Hamiltonian describing the TJJ [22], which we now 
briefly describe in terms of its bosonic and fermionic de- 
grees of freedom. 

Description of the junction. — The following Hamilto- 
nian governs the dynamics of the relative phase degree of 
freedom across the junction [32] 

(1) 

where ip(x) is the relative phase, x is the coordinate run- 




FIG. 1. (Color online) An annular topological Josephson junc- 
tion trapping a single soliton (green). Counter-propagating 
Majorana edge states (red and blue lines) are nucleated at 
the junction. A charge Q and phase $ are induced externally 
within the central region. 
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ning along the junction, A the Josephson penetration 
length, c the renormalized velocity of light, and g a di- 
mensionless constant which depends on the parameters 
of the junction. In addition, the junction supports two 
counter-propagating Major ana edge states with a Hamil- 
tonian = H + H (H describes the external edge and 
H the internal) 



H = i- I dxi/j(x)d x i/j(x), 



H 



=-'iJ 



dxifj(x)d x ilj(x). 



(2) 



Here x G [0, L] runs around the circumference of the edge 
and v is the neutral edge velocity. In addition, there ex- 
ists a tunneling term between the two edge states which 
couples the phase degree of freedom and the two Majo- 
rana degrees of freedom, taking the form 



H t/I 



dx W \x) / i/j(x)i/j(x) , 



(3) 



where W(x) = mcos [cp(x)/2] is the Majorana mass term 
[2Q| [22] . The full Hamiltonian for the junction, + 
H^p + H tU n [H], is an extension of the supersymmetric 
Sine-Gordon model for general values of m [26 . 
We perform the following mode expansion, 



v>0) = J j J2 



-2ninx/ L 



4>n, 



(4) 

The fields obey anti-commutation relations of the form 
Mx), (*')} = ± £ e~ 2 ^-^l L = 5(x - x'), 

n 

where the modes ip n satisfy {V^Vv} = <Wn',o (and 
similarly for the opposite chirality). In addition, 
{ip(x), ^(x f )} = 0. Note that in particular this implies 
2ipQ = 1, or ipQ = 1/2 (for either chirality). Plugging this 
into the Hamiltonian we get, 



H 



H 



2ttv 
2ttv 



n 



2ttu 
L ^ 



2irv 



C. 



(5) 



It is the quantity C (and £), the dimensionless momen- 
tum operator, that we will pursue in the following. Using 
Eq. Q, periodic boundary conditions on the Majorana 
field imply n G Z, while anti-periodic boundary condi- 
tions imply n G Z + \ . 

A quantized shift of the edge states ' total momentum. — 
To understand the change in the momentum when the 



boundary conditions are exchanged between periodic and 
anti-periodic, we write C as 



n>0 



2 ^ 

n>0 



(6) 



n>0 



where Co is a constant, which can be thought of as the 
ground state contribution; the other part consists of ex- 
citations above the ground state. If there is a vortex 
enclosed by the edge, n G Z, otherwise n G Z + 1/2. 

For a closed circular Josephson junction, in the absence 
of tunneling, we have the following 



C = ^2 n^-n^n + £ (A^mod2), 



n>0 

£ 

m>0 



rmp-mipm + C (N v mod2), 



(7) 
(8) 



where N v (N v ) denotes the number of vortices enclosed 
by the external (internal) edge. We now calculate the 
difference in the ground state contribution in the presence 
of a Josephson vortex within the junction, i.e. N v = 1 
and N v = 0. We employ a regularizing function F(x) 
such that F'(x) = d x F(x) decays to zero faster than 
1/x 2 when x — » oo, and F'(0) = 1. We calculate the 
regularized sum [33] 



AC = £ (1) - A)(0) 



(9) 



2 ^ 

n=l 



nF '(an) — (n 
By taking the limit a — >• we now get 



l -)F'{a{n-\)) 



F{an) - F(a(n - -)) 



1 CO 

Z n=l 

i 00 r 2i 



n=l 
poo 



d a / d(an) 



a/2 



a 



~-F'(ari) - -F"{an) 



^o )+ nx>)] = i. 



(10) 



This result is well-known from conformal field theory (see 
e.g. [34]) where a zeta function regularization is typically 
used to obtain the same value. In the following we ex- 
plore how this quantized shift can affect the dynamics of 
the soliton, which is governed by field theories that are 
explicitly non-conformal. 

Effective Hamiltonian for a Josephson vortex. — We 
consider the solution for a classical soliton (in the non- 
relativistic limit and for A<L) 



(p 8 (x, q(t)) = 4 arctan 



exp 



(ii) 
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with a center of mass coordinate at q{t). We plug the so- 
lution into the Lagrangian, using d t (p s — QdxPs, and the 
integrals dx (d x (p s ) 2 = | and dx (1 — cos<^ s ) = 
4A to get 



:m s q 



(12) 



where we defined the soliton mass m s and the soliton 
rest energy. We now proceed to the Major ana sector, 
with ^ = ( ij) 



where 



Hn 



ivd x iW{x,q{t)) 
-iW{x,q{t)) —ivd x 



(13) 



(14) 



and 



W(x,q(t)) 



mcos[tp s (x,q(t))/2] 
mtanh [{x — q{t)) /A] 



W(x 



(15) 



The equations simplify considerably by taking a boost to 
the moving frame 



x = x 
t' t, 
d x = dx' 
d t = ~qd x 



(16) 



We see that the Major ana fields couple to the center of 
mass velocity of the soliton via a vector-potential like 
term that measures the total momentum carried by the 
two counter-propagating edge states, taking the form: 

/_ _ 2tt 
dx {^d x ^j + $d x $) = —q(£ - C). (17) 

The full Hamiltonian + + H tun , written in the 
background of a single soliton, is given in terms of the 
soliton's center of mass momentum p (which we now re- 
instate as a quantum operator) as 

l 2 



Ha — En 



2m. 



P 



2tt 
~L 



(£-£) 



C) + i / dxW(x)ip(x)<i/j(x). (18) 



2ttv , „ 

+ — ( £ 

This Hamiltonian describes the dynamics of the Joseph- 
son vortex within the junction and is a central result of 
this paper. Not only it captures universal features as 
we now describe, it can also be used to extract potential 
deviations from universality. 

The Hamiltonian admits a subtle symmetry, that the 
edge states' total momentum operator and tunneling op- 
erator, given respectively by 

2-7T 

P = -(£-£), (19) 



dx W(x)ip(x)i/;(x), 



(20) 



in fact commute. To show this, one can use the following 
argument. First, we calculate the commutator 



[P,T] 



2tt 

T 



dx d x W(x)ip(x)ijj(x). 



(21) 



For our choice of </? s , the term d x W(x) is localized around 
the center of the soliton. We can formally write the com- 
mutator in Eq. (1211) as [P,T] =J^\im x ^ Xo d Xo T[x ], 



where T[xq] is obtained from Eq. (20) by replacing W(x) 
by W(x — xn). Since T[xn) does not depend on xn owing 
to translational symmetry, the derivative term is zero. 
Hence [P, T] = 0. The free spectrum of the operator P, 
given by 



2tt 

T 



1 

16 



n e Z, 



(22) 



remains unaffected by the tunneling, with ± related to 
the parity of the number of vortices within the path. This 
demonstrates a certain robustness of the 1/16 result even 
for this non-conformal case. 

Proposed experiment for detecting the phase shift. — 
In the absence of a Josephson vortex within the junc- 
tion, the tunneling results in the Majorana modes gap- 
ping out. When a Josephson vortex is then introduced, a 
zero energy bound state nucleates on the vortex while the 
rest of the modes remain gapped. There are two cases 
of interest: (i) If there is an even number of vortices in 
the central region of the sample, the total dimension- 
less momentum carried by the two counter-propagating 
edge states is 1/16; and (ii) If there is an odd number 
of vortices the same momentum is —1/16. Hence, while 
the total energy for the two edges is the same for the two 
cases, a certain part of the momentum operator, £0 — A), 
changes sign. 

The energy spectrum of the Josephson vortex can be 



derived from Eq. (18), and in the presence of an exter- 



nally induced Aharonov-Casher charge Q within the cen- 
tral region, is given by 



Eft — En 



Q 

2e 



(V± + !M 

V 4 16/ 



n b 



(23) 



where E c is the charging energy, rif = ± is the fermion 
parity within the enclosed path of the Josephson vortex, 
n v = ± the parity of the number of vortices within the 
same region, and G Z the relative number of Cooper 
pairs between the two superconducting plates. In the low 
energy sector there is an emergent dependence between 
nf and n v : If n v = 1 then nf = 1, but if n v = —1 then 
rif is free [22] . 

We start from the case rif = 1 and n v = 1, tuning the 
junction into its zero resistance state by shifting the in- 
duced Aharonov-Casher charge Q; the Josephson vortex 
accordingly acquires a vanishing velocity. Next we add 
an extra vortex within the central region of the sample, 
shifting the value of n v to —1. The Josephson vortex then 
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FIG. 2. (Color online) Energy spectrum for the Josephson 
vortex. Dashed green lines describe the energy of the Joseph- 
son vortex in the presence of an even number of vortices 
enclosed within its path. Solid blue and red lines describe 
the case with an odd number of vortices (for even and odd 
fermion parities respectively). The shift between these two 
cases [dashed vs solid] is a universal 1/4 (in units of the elec- 
tronic charge e). 

performs a persistent motion and the junction is driven 
into its finite resistance state (as measured between the 
internal and external superconductor plates). As men- 
tioned before, we now have two knobs to control the per- 
sistent motion of the vortex: (i) The induced charge Q 
(through the Aharonov-Casher effect); and (ii) the flux 
<£, quantized in units of h/2e. Each additional vortex 
leads to a phase shift which is equivalent to a ±e/4 shift 
in the induced Aharonov-Casher charge (see Fig. |2| 

Discussion. — A central result of this paper is the 
identification of a relative 7r/4 phase associated with a 
Josephson vortex in a topological Josephson junction en- 
circling an odd versus even number of vortices. It is use- 
ful to compare this result with the full conformal case 
captured by a standard fusion rule from conformal field 
theory (see, e.g., [34), a(z)a(0) ~ z' 1 ^ 8 [I + z l / 2 ^(z)] . 
The last fusion rule is associated with a single chiral Ma- 
jorana edge state, with / the identity field and xp and 
a fields of dimensions 1/2 and 1/16 respectively. If we 
identify the field a{z) as the vortex and z = x + iy as its 
coordinate, then this equation reproduces the presence 
of a — 7r/4 phase shift for a rotation of one vortex around 
another, z — >• e 27TZ z. For the case of an odd fermionic 
number, a 37r/4 phase shift would ensue. Our theory is 
explicitly non-conformal, yet the topological data decided 
by the conformal case is fully reproduced here, demon- 
strating the rigidity of the underlying topological quan- 
tum field theory. Finally, we address the potential effect 
of several experimental factors that were not explicitly 
accounted for in the theory described above. First, inter- 
actions between the vortices are not expected to produce 
geometric phases; other effects of interaction are deferred 



to a later work. Second, quasi-particle tunneling may af- 
fect the part of the spectrum that is sensitive to fermion 
parity effects, however, the exact e/4 shift discussed here 
remains immune to a shift by e, hence so is the residual 
motion of the soliton generated by it. 
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